Origin and roles of a strong electron-phonon interaction 
in cuprate oxide superconductors 
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A strong electron-phonon interaction arises from the modulation of the superexchange interaction 
by lattice vibrations. It is responsible for the softening of the half-breathing modes around (±7r/a, 0) 
and (0, ±7r/a) in the two-dimensional Brillouin zone, with a being the lattice constant of CUO2 
planes, as is studied in Phys. Rev. B 70, 184514 (2004). Provided that antiferromagnetic spin 
fluctuations are developed around Q = (±37r/4a, ±7r/o) and (±7r/a, ±37r/4a), the electron-phonon 
interaction can also cause the softening of Cu-O bond stretching modes around 2Q, or around 
(±7r/2a, 0) and (0, ±7r/2a). The softening around 2Q is accompanied by the development of charge 
fluctuations corresponding to the so called 4a-period stripe or 4a x 4a-period checker-board state. 
However, an observation that the 4a-period modulating part or the 2Q part of the density of states is 
almost symmetric with respect to the chemical potential contradicts a scenario that the stabilization 
of a single-2Q or double-2Q charge density wave following the complete softening of the 2Q bond 
stretching modes is responsible for the ordered stripe or checker-board state. It is proposed that 
the stripe or checker-board state is simply a single-Q or double-Q spin density wave, whose second- 
harmonic effects can explain the observed almost symmetric 2Q part of the density of states. The 
strong electron-phonon interaction can play no or only a minor role in the occurrence of d7-wave 
superconductivity in cuprate oxides. 

PACS numbers: 71.38.-k, 74.20.-z, 75.30.Et 



I. INTRODUCTION 

It is one of the most interesting and important is- 
sues in condensed-matter physics to elucidate the mech- 
anism of high critical temperature (high-T c ) supercon- 
ductivity occurring in cuprate oxides>i The oxides are 
highly anisotropic quasi-two-dimensional oxides, whose 
main compositions are Cu02 planes. High-T c super- 
conductivity occurs on the CuC>2 planes. There are 
pieces of evidence that the electron-phonon interaction 
is strong on the CuC>2 planes: the softening of the half- 
breathing modes around (±7r/a, 0) and (0, ±7r/o) in the 
two-dimensional Brillouin zone (2DBZ)»2iMi5i& with a be- 
ing the lattice constant of the CuC>2 planes, the softening 
of Cu-0 bond stretching modes around (±7r/2a, 0) and 
(0, ±7r/2a) in 2DBZ,- 7 ' 8 kinks in the dispersion relation of 
quasiparticles j 9 -^ and so on. It may be argued therefore 
that the electron-phonon interaction must play a major 
role in the occurrence of high-T c superconductivity. On 
the other hand, observed isotope shifts of T c are smalljii 
which implies that the strong electron-phonon interaction 
can play only a minor role in high-T c superconductivity 
itself. The origin and roles of the strong electron-phonon 
interaction should be clarified in order that the issue of 
high-T c superconductivity might be solved. 

Parent cuprate oxides with no doping are Mott in- 
sulators. When holes or electrons are doped into the 
the Mott insulators, high-T c superconductivity appears. 
Cuprate oxide superconductors lie in the vicinity of the 
Mott metal- insulator transition or crossover. It may be 
argued therefore that strong electron correlations must 
play a crucial role not only in the occurrence of high-T c 
superconductivity but also in the origin and roles of the 



strong electron-phonon interaction. 

The Hubbard model is one of the simplest effective 
Hamiltonians for strongly correlated electron liquids. In 
Hubbard's approximation ; 12 ! 13 a band splits into two 
subbands when the on-site repulsion U is so large that 
U > W , with W being the bandwidth of unrenormalized 
electrons. The subbands are called the upper Hubbard 
band (UHB) and the lower Hubbard band (LHB), and a 
gap between UHB and LHB is called the Hubbard gap. 
In Gutzwiller's approximation ) 1 ! 15 ' 16 a narrow quasipar- 
ticle band appears around the chemical potential. The 
band and quasiparticles are called the Gutzwiller band 
and quasiparticles, respectively. It is plausible to spec- 
ulate that the density of states has in fact a three- 
peak structure, with the Gutzwiller band between UHB 
and LHB. Both of the approximations are single-site ap- 
proximations (SSA). Another SSA theory confirms the 
speculation t — showing that the Gutzwiller band appears 
at the top of LHB when the electron density per unit cell 
is less than one. The nature of the ground state of the 
Hubbard model depends on the nature of the Gutzwiller 
quasiparticles. 

The SSA that considers all the single-site terms is re- 
duce to determining and solving self-consistently the An- 
derson model ) 18 ' 19 ' 20 which is one of the simplest effec- 
tive Hamiltonian for the Kondo problem. Hence, the 
three-peak structure corresponds to the Kondo peak be- 
tween two subpeaks in the Anderson model, or in the 
Kondo problem. The s-d model is also one of the sim- 
plest effective Hamiltonian for the Kondo problem. Ac- 
cording to Yosida's perturbation theory 2 ^ and Wilson's 
renormalization-group theory^ the ground state of the 
s-d model is a singlet or a normal Fermi liquid provided 
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that the Fermi surface of conduction electrons is present. 
Since the s-d model is derived from the Anderson model, 
the ground state of the Anderson model is also a normal 
Fermi liquid. It is certain therefore that under the SSA 
the ground state of the Hubbard model is a normal Fermi 
liquid or a metali^ 3 . Even if the Hubbard gap opens, the 
Fermi surface of the Gutzwiller quasiparticles is present. 

The SSA can also be formulated as the dynamical 
mean-field theory 2 ^ and the dynamical coherent poten- 
tial approximation^ In the SSA, local fluctuations are 
rigorously considered but Weiss mean fields, which are 
responsible for the appearance of the corresponding or- 
der parameter, are ignored. Hence, the SSA is rigorous 
for infinite dimensions within the Hilbert subspace with 
no order parameter^ In Kondo-lattice theor y 18 ' 19 ' 20 an 
unperturbed state is the normal Fermi liquid, which is 
constructed in the non-perturbative SSA theory, and not 
only effects of intersite fluctuations but also ordering due 
to Weiss mean fields such as magnetism, superconduc- 
tivity and so on are perturbatively considered. Kondo- 
lattice theory can also be formulated as 1/d expansion 
theory, with d being the spatial dimensionality. 

The d-p model, where Cu 3c? and O 2p orbits are ex- 
plicitly considered, is one of the simplest effective Hamil- 
tonians for cuprate oxide superconductors. Since the on- 
site repulsion U plays a crucial role in the d-p model as 
it does in the Hubbard model, it is straightforward to 
extend the analysis for the Hubbard model to the d-p 
model. Observed quasiparticle states, which are often 
called mid-gap states, are simply the Gutzwiller quasi- 
particle states, which can also be renormalized by inter- 
site fluctuations. When observed specific heat coefficients 
as large as 7 ~ 14 mJ/K 2 mol are used) 27 ! 28 the Fermi- 
liquid relation give s 29 ! 30 

W* — 0.3-0.4 eV, (1.1) 

for the effective bandwidth of the Gutzwiller quasipar- 
ticles in optimal-doped cuprate oxide superconductors, 
where T c is the highest as a function of doping concen- 
trations. According to field theory, the superexchangc 
interaction arises from the virtual exchange of a pair 
excitation of electrons across the Hubbard gapj^I Since 
the Gutzwiller quasiparticles, which are responsible for 
metallic properties, plays no significant role in the virtual 
exchange process, the superexchange interaction is rele- 
vant even in a metallic phase, provided that the Hubbard 
gap opens. Cooper-pairs can also be bound by a mag- 
netic exchange interaction^ Since the superexchange in- 
teraction constant is as strong as^ 3 . J = —(0.10-0.15) eV 
between nearest-neighbor Cu ions on a Cu02 plane, ob- 
served high T c can be easily reproduced. In actual, it 
has already been proposed that the condensation of dj- 
wave Cooper pairs between the Gutzwiller quasiparticles 
due to the superexchange interaction is responsible for 
high-T c superconductivity ! 34 ' 35 Since the superexchangc 
interaction is strong only between nearest-neighbor Cu 
ions, it is definite that theoretical T c of the d'y wave is 
much higher than those of other waves. In fact, high- 



T c superconductivity occurs in an intermediate-coupling 
regime | J\/W* = 0.3-0.5 for superconductivity, which is 
realized in the strong-coupling regime for electron corre- 
lations defined by U/W > 1. 

Since charge fluctuations are suppressed by strong elec- 
tron correlations, the conventional electron-phonon inter- 
action arising from charge-channel interactions must be 
weak in cuprate oxide superconductors. On the other 
hand, an electron-phonon interaction arising from spin- 
channel interactions can be strong. For example, an 
electron-phonon interaction arising from the modulation 
of a magnetic exchange interaction by phonons plays a 
significant role in the spin-Peierls effect. It has been 
shown in a previous paper— that an electron-phonon 
interaction arising from the modulation of the superex- 
change interaction by phonons is strong in cuprate oxide 
superconductors. The electron-phonon interaction can 
explain the softening of the half-breathing modes around 
(±7r/o, 0) and (0,±7r/a) in 2DBZ. It has been predicted 
that the softening must be small around (±7r/a, ±7r/o) 
in 2DBZ. An attractive mutual interaction due to such 
an electron-phonon interaction is strong between quasi- 
particles on next-nearest-neighbor Cu ions, but is very 
weak between those on nearest-neighbor Cu ions. There- 
fore, the mutual interaction can play no significant role 
in the binding of cfy-wave Cooper pairs. Observed small 
isotope shifts of T c can never contradict the presence of 
the strong electron-phonon interaction. 

The so called 4a-period stripes and 4a x 4a-period 
checker boards are observed in under-doped cuprate ox- 
ide superconductors ) 37 ' 38 ! 39 ' 40 ! 41 ' 42 whose doping concen- 
trations are smaller than those of optimal-doped ones. 
The wave numbers of Cu-0 bond stretching modes, 
(±n/2a, 0) and (0, ±n/2a) in 2DBZ, correspond to the 
period 4a of stripes and checker boards. The softening 
of the stretching modes is accompanied by the develop- 
ment of 4a-period or 4ax4a-period fluctuations in charge 
channels, which are simply stripe or checker-board fluc- 
tuations. It may be argued therefore that a charge den- 
sity wave (CDW) following the complete softening of the 
bond stretching modes is responsible for ordered stripes 
and checker boards. 

One of the purposes of this paper is to show that the 
strong electron-phonon interaction can also explain the 
softening of Cu-0 bond stretching modes in cuprate ox- 
ide superconductors. The other purpose is to examine 
critically the relevance of the CDW scenario, whether or 
not the CDW is actually responsible for ordered stripes 
and checker boards. This paper is organized as fol- 
lows: Preliminary discussions are presented in Sec. [TTJ 
the derivation of the electron-phonon interaction is re- 
viewed in Sec. Ill Al and Kondo-lattice theory is reviewed 
in Sec. Ill Bl The softening of the bond stretching modes 
around (±7r/2a, 0) and (0,±7r/2a) in 2DBZ is studied in 
Sec. IIIII The relevance of the CDW scenario for stripes 
and checker boards is critically examined in Sec. II VI An 
argument on the mechanism of high-T c superconductivity 
is given in Sec. [V] Conclusion is presented in Sec. IVII 
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II. PRELIMINARIES 

A. Electron-phonon interaction 

In cuprate oxide superconductors, the superexchange 
interaction arises from the virtual exchange of a pair ex- 
citation of 3d electrons between UHB and LHB that are 
strongly hybridized subbands between Cu 3d and O 2p 
orbits^i. When the broadening or finite bandwidths of 
UHB and LHB are ignored, the exchange interaction con- 
stant between nearest-neighbor Cu ions on a CuC>2 plane 
is given by 



J = 



4U 4 



1 



1 



(e d -e p + U) 2 \e d -e p + U U J ' 



(2.1) 



with V being the hybridization matrix between nearest- 
neighbor O 2p and Cu 3d orbits, e d and e p the depths 
of Cu 3d and O 2p levels, and U the on-site repulsion 
between Cu 3d electrons. 

Doped holes reside mainly at O ions. The preferential 
doping suggests that O 2p levels are shallower than Cu 3d 
levels or that parent cuprate oxides with no hole doping 
must be charge-transfer insulators rather than Mott in- 
sulators; charge-transfer insulators and Mott insulators 
are characterized by e p > e d and e p < e d , respectively. 
Since the hybridization between Cu 3d and O 2p orbits 
must be strong, it may also be argued that Cu 3d lev- 
els are much deeper than O 2p levels, that is, e p ^> e d 
rather than e p > e d , to explain the observed preferential 
doping. However, the suggested level scheme of e p > e d 
or e p ^> e d disagrees with the prediction of Mott insu- 
lators, e d — e p ~ 1 eV, by band calculations i 43 i 44 ' 45 The 
preferential doping does not necessarily mean that the 
parent cuprate oxides are charge-transfer insulators, but 
it simply means that the local charge susceptibility of 
3d electrons is much smaller than that of 2p electrons, 
which implies that the effective on-site repulsion U be- 
tween 3d electrons is very strong. It is assumed in this 
paper that V ~ 1.6 eV and e d — e p ~1 eV, as is pre- 
dicted by band calculationsj 43 ' 44 i 45 Since the on-site U 
should be so large that the Hubbard gap might open, 
it is assumed that U ~ 5 eV. Then, Eq. (|2.1[) gives 
J ~ —0.27 eV. This is about twice as large as the ex- 
perimental J = —(0.10-0.15) eVi 3 - 3 - This discrepancy is 
resolved when nonzero bandwidths of UHB and LHB are 
considered^ 1 - 

Displacements of the ith Cu ion and the [ij]th O ion, 
which lies between the nearest-neighbor ith and jth Cu 
ions, are given by 



and 



j£ y/2NM d Lo Xci 



j£ ^2NM p u Xc 



- Aq 



(2.2) 



(2.3) 



with Hi and R[ij] = (R; + Rj)/2 being positions of the 
ith Cu and [ij]th O ions, M d and M p masses of Cu and O 
ions, b\ and 6> q creation and annihilation operators of 
a phonon with a polarization A and a wave vector q, LO X q 
a phonon energy, eA q = (eAq £Aq,y, £Aq, 2) & polarization 
vector, and N the number of unit cells. Here, only lon- 
gitudinal phonons are considered so that it is assumed 
that 6A q = (fe, %, <Zz)/|q| for q within the first Brillouin 
zone. The q dependence of v d . x<i and w P ,Aq is crucial. 
For example, t>d,Aq = and v p ^ X q = 0(1) for modes that 
bring no change in adjacent Cu-Cu distances. 

Two types of electron-phonon interactions arise from 
the modulations of the superexchange interaction J by 
the vibrations of O and Cu ions. When they are con- 
sidered, it is convenient to define a dual-spin operator. 
First, a single-spin operator is defined by 



S(q) 



—y 1 - 



°' a ^ d (k+iq)a d (k-iq)/3' 



(2.4) 



with er = (<r x , a y , a z ) being the Pauli matrixes and d^ ka 
and dko- being creation and annihilation operators, re- 
spectively, of 3d electrons with wave number k. Then, 
the dual-spin operator is defined by 

^r(q) = i$>r(q')[S(q'+iq)-S(-q'+i q )] , (2 . 5 ) 



q' 



with 



and 



r] s (q) = cos(q x a) + cos(q y a), 



Vd(q) = cos(q x a) - cos(q y a). 



(2.6) 



(2.7) 



It is assumed in this paper that the x and y axes are 
within C11O2 planes and the z axis is perpendicular to 
CuC>2 planes. The electron-phonon interactions are sim- 
ply given by 



n P = ic p ^ 



^ y/2NMpUJXq 

x%(q) Vr (|q)Pr(q), 



(2.8) 



r=s.d 



and 



H d = iC d Y, 



y/2NM d U X q 



^-q + ^Aq 



*?r(q)Pr(q), 



(2.9) 



r=s,d 



with C p and C d being real constants, which are given in 
the previous paper^ and 



■n s (q) = 2 \e Xx sin I — J + e Xy sin J 
%(q) = 2 e Xx sin (ppj ~ e Xy sin (^pj 



(2.10) 
(2.11) 
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B. Kondo-lattice theory 

One of the simplest effective Hamiltonians for the elec- 
tron part of cuprate oxide superconductors is the d-p 
model on a square lattice. Since the anisotropy is large, it 
is convenient to consider phenomenologically quasi-two- 
dimensional features. The d-p model is approximately 
mapped to the t-J model or the t- J-infinite-?7 model4£ 



H 



t-J = ^Ujd^dja - ~J ^(S, • Sj) 



+Uoo ^ (f.^difdl^dii 



(2.12) 



with the summation over (ij) indicating that the sum- 
mation should be made over nearest neighbors and 



s 4 = ]T^dU 

a/3 



i/3- 



(2.13) 



The carrier density per unit cell is defined by 

^ (2.14) 



It should be noted that the infinitely large on-site repul- 
sion Uoo is introduced to exclude double occupancy so 
that n can never be larger than unity, or < n < 1. 
The electron and hole pictures should be taken for the 
so called hole-doped (n < 1) and electron-doped (n > 1) 
cuprate oxide superconductors, respectively, so that n is 
the electron density for hole-doped ones and is the hole 
density for electron-doped ones. The doping concentra- 
tion is defined by 6 = 1 — n, and the optimal concentra- 
tion, where superconducting T c is the highest as a func- 
tion of S, is S ~ 0.15. Then, 5 < 0.15 and S > 0.15 are 
called under-doped and over-doped concentrations, re- 
spectively. When transfer integrals between nearest and 
next-nearest neighbors, which are denoted by t and t' , 
are only considered, the dispersion relation of electrons 
or holes is given by 



£(k) 



+4t' cos(k x a) 



2t [cos(k x a) + cos(kyd)] 
) cos(k y a) 



(2.15) 



According to band calculations ; 43 ' 44 ' 45 it follows that t = 
— (0.3-0.5) eV and t' ~ — 0.3t for electrons in hole-doped 
cuprate oxide superconductors and t = —(0.3-0.5) eV and 
t' ~ +0.3t for holes in electron-doped ones. 

Every physical quantity is divided into single-site and 
multi-site terms. Calculating the single-site term is re- 
duced to determining and solving selfconsistently the An- 
derson model, as is discussed in the Introduction. When 
it is assumed that there is no order parameter, for exam- 
ple, the self-energy of electrons is divided into single-site 
and multi-site self-energies: 



£ CT (i£„,k) = £ CT (ie„) + AS (T (ie„,k). 



(2.16) 



The single-site self-energy S CT (ie n ) is given by that of the 
Anderson model. It is expanded as 



1- 



+ (l-4> s )-crg f i B H + 0(e 2 ), (2.17) 

at T = K in the presence of an infinitesimally small 
Zeeman energy g/isH, with g the g factor and /jlb the 
Bohr magneton. The expansion coefficients S , </>_,, and 
4> s are all real; (f> 8 ~ 2</> 7 > 1 for n ~ 1. When Eq. (|2.17|) 
is used and the multi-site self-energy is ignored, the dis- 
persion relation of the Gutzwillcr quasiparticles is given 
by 



UK = t 



So + E(k) - (i 



^aW s g^ B H, (2.18) 



with 



W s = 4> s /4>-f, 



(2.19) 



being the so called Wilson ratio for the Kondo problem. 

The irreducible polarization function in spin channels 
is also divided into single-site and multi-site polarization 
functions: 



tt s (iui , q) = 7r s (iuji ) + Att s (iu>i , q) . 



(2.20) 



The single-site polarization function n s (iuii) is given by 
that of the Anderson model. The spin susceptibilities of 
the Anderson and t-J models are given, respectively, by 



and 



1 - UooTTs^ilVl) ' 



with 



J(q) = 2 J [cos(g x o) 



<s(q y a)} 



(2.21) 



(2.22) 



(2.23) 



In Eqs. (|2.2ip and (|2.22p . the conventional factor \g 2 p? B 
is not included. A physical picture for Kondo lattices 
is that local spin fluctuations at different sites interact 
with each other by an intersite exchange interaction. In 
Kondo-lattice theory, therefore, an intersite exchange in- 
teraction I s (iu>i , q) is defined by 



X s (*w ; ,q) 



It follows that 



I s (iui,q)x s (iuJi) 



I s (iwi,q) = J(q) +2£/^A7r s (iw ; ,q). 



(2.24) 



(2.25) 



The derivation of Eq. ([2~25]) from Eqs. (|2~2T|) and 
(I2.22[) is rigorous because ignored terms, which are 
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0[l/UooXs{iwi)], vanish for infinitely large Uoo- The term 
of 2U^ ATT s (iu)i : q) is composed of two terms4£ 

2f/^A7r s (iwi,q) = J Q (m, q) - 4A(iw,, q). (2.26) 

The first term Jq (icj; , q) is an exchange interaction aris- 
ing from the virtual exchange of a pair excitation of 
the Gutzwiller quasiparticles. According to the Ward 
relation^ the static component of the single-site irre- 
ducible vertex function in spin channels is given by 

\ s =&[1-E7oo7r.(0)]. (2.27) 
Then, it follows that 



UooXs = 24> s /xs(0). 



(2.28) 



When the vertex correction X s given by Eq. (|2.28[) is used, 
it follows that 



J Q {iu h q) 
with 

with 



N ^ 



/K g (k)]-/[g ff (k + q)] 
C CT (k + q)-^(k)+^; ! 



/(£) = 



1 



(2.29) 



(2.30) 



(2.31) 



gE/fcfiT I J ' 

In Eq. (|2.29[) . the single-site term is subtracted because 
it is considered in the SSA. The strength of the ex- 
change interaction is proportional to the width of the 
Gutzwiller band. Since the chemical potential lies around 
the center of the Gutzwiller band or the nesting of the 
Fermi surface is sharp, the exchange interaction is an- 
tiferromagnetic (AF) in cuprate oxide superconductors. 
The second term — 4A(jw/,q) corresponds to the mode- 
mode coupling term in the self-consistent renormaliza- 
tion theory of spin fluctuations^ which is relevant in the 
weak-coupling regime for electron correlations defined by 
U/W < 1. 

When Eq. (|2 . 28[) is used, the mutual interaction medi- 
ated by intersite spin fluctuations, which works between 
the Gutzwiller quasiparticles, is given by 

(£/ooA s ) 2 [ Xs ^,q) -xs(m)] = 4> 2 s I* s (iu> h q), (2.32) 



with 



^ (m, q) 



J s (iwz,q) 



1 - jl s {iuj u c0x s {i^i) 



(2.33) 



In Eq. (|2.32p . the single-site term is subtracted because 
it is considered in the SSA. The exchange interaction 
I a (iuii, q) is enhanced into I*(iu)i, q) by spin fluctuations. 
The mutual interaction mediated by spin fluctuations is 
essentially the same as that due to the exchange interac- 
tion I s (iuji,q) or I*(iui,q). 

In Kondo-lattice theory, a unperturbed state is con- 
structed in the non-perturbative SSA theory and multi- 
site or intersite effects are perturbatively considered in 
terms of I s (iu>i, q) or Ig(iu>i, q). 



III. SOFTENING OF PHONONS DUE TO 
ANTIFERROMAGNETIC SPIN FLUCTUATIONS 

An effective Hamiltonian to be eventually examined in 
this paper is 



with 



H = Tit- j + Ttph + Tt p + Hd, (3.1) 



Wph = £ W Aq( 6 W*q+£), (3-2) 
Aq V ' 



and Hp and Hd defined by Eqs. (|2.8[) and (|2.9[) . respec- 
tively. The t-J model Ht-j is defined on a square lat- 
tice in Sec. Ill Bl but it should be defined on a quasi-two- 
dimensional lattice here. Although no interlayer coupling 
is included in Ht-j, the nature of quasi- two-dimensional 
AF spin fluctuations is phenomenologically considered, 
which plays one of the most crucial roles in the softening 
of Cu-0 bond stretching modes, as is examined below. 

In the absence of H v and Hd-, the Green function for 
phonons is given by 



2uj 



Aq 



(icvi) 2 - bJ 2 



(3.3) 



Aq 



When the self-energy for phonons is denoted by 

^a ' l) ' wn i cn can be perturbatively calculated in 
terms of H p and Hd, the renormalized Green function for 
phonons is given by 

D x (iu)i,q) = D^'(iu; h q) 



-D { °\iw h q)Z^ h) (m,q)D x (iu l ,q) 



J A 
2WAq 



{iuji) 2 - uj 2 ^ - 2u Aq £ A ph) (iwz, q) 



(3.4) 



The renormalized energy of phonons, which is denoted 

by w Aq > is s iven b y 

K q ) 2 - ^ 2 q - 2^ Aq 4 ph) ( W * q + iO, q) = 0. (3.5) 
The renormalization of phonon energies is given by 



Aw Ac 



= E 



; Aq _ ^Aq 
A ph) (c Aq + zO,q) 



[4 ph) ( WAq + i 0,q)] Z /2^ q + ... , (3.6) 



unless |S A ph ^(w Aq + iO, q)| is larger than |w Aq |. 

As well as AF spin fluctuations are developed around 
Q due to I s (iu>i,q), dj-wave superconducting (SC) and 
charge bond order (CBO) fluctuations are also developed 
due to I s (iu>i,q) or I*(iuJi,q). Although charge fluctua- 
tions are never much developed, charge-channel fluctu- 
ations can also contribute to the softening of phonons, 
as well as AF, SC, and CBO fluctuations, provided that 
vertex corrections for the dual spin operator in spin, SC, 
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and CBO channels are properly treated. According to 
the previous paper,— the softening of the half-breathing 
modes is mainly caused by the charge-channel fluctua- 
tions. Since the charge-channel fluctuations are signifi- 
cant in the metallic phase, the softening is large in the 
metallic phase but is small in the insulating phase. 

Since phonons can couple with two lines or two chan- 
nels of spin fluctuations to the lowest or first order in 
the dual-spin operator, as is shown in Eqs. (|2.8[) and 
(|2.9p . AF spin fluctuations around (±37r/4a, ±7r/a) and 
(±7r/a, ±37r/4a) in 2DBZ can play a significant role in 
the softening of Cu-0 bond stretching modes around 
(±7r/2a,0) and (0,±7r/2a) in 2DBZ. According to Mer- 
min and Wagner^ if the Neel temperature Tjv were 
nonzero in two dimensions integrated effects of two- 
dimensional critical AF spin fluctuations would be di- 
vergent at TV, which leads to a conclusion that T/v 
must be zero in two dimensions. Their argument im- 
plies that quasi-two-dimensional critical AF spin fluctu- 
ations can play a crucial role in the softening, at least, 
in an AF critical region of cuprate oxide superconductors 
provided that the anisotropy of critical AF spin fluctua- 
tions is large. In order to examine how crucial a role the 
anisotropy plays in the softening, it is more convenient to 
use a phenomenological expression for the spin suscepti- 
bility, which includes explicitly the anisotropy factor for 
AF spin fluctuations, than to calculate microscopically 
the spin susceptibility for quasi-two-dimensional systems. 
The superexchange interaction J(q), which is given by 
Eq. (|2~33"|) . has broad peaks at (±n/a,±n/a) in 2DBZ 
and the exchange interaction Jq(0, q), which is given by 
Eq. (|2.29[) . has sharp peaks at nesting wave numbers of 
the Fermi surface. Therefore, it is assumed in this paper 
that J s (0,q) is maximal at Q = (±37r/4a, ±7r/a, Q z ) and 
(±n/a, ±37r/4a, Q z ) or that x s (0, q) is maximal at Q, 
and that the spin susceptibility (|2.24|) is approximately 
but well described by 



Xs (i^,Q + q) 



Xs(0,Q)k 2 



K 2 + ( qil a) 2 + 5 2 (q z c) 2 + 



Faf 



(3.7) 



around each of Q's, with q|| = (q x ,q y ) being the com- 
ponent parallel to CuC>2 planes, q z the component per- 
pendicular to CuC>2 planes, c the lattice constant along 
the z axis, and Taf the energy scale of AF spin fluctua- 
tions. The anisotropy factor 5 is introduced to consider 
quasi-two-dimensional AF spin fluctuations. The corre- 
lation length within the x-y plane is ajn and that along 
the z axis is 8c/ k. A cut-off q c = ir/3a is introduced 
in such a way that \s (i^>i, Q + q) = for \q x \ > q c or 
1^1 > q c . The anisotropy of the lattice constants plays 
no role when 6 and q c arc defined in these ways. 

When AF spin fluctuations are only considered, the 
self-energy for phonons is given by 3 - 6 - 



4 ph Wq) 



2M p w Aq . |r 
xA rr '(iw;,q) 



-5>(q)r r ,(q) 



with 

>r(q) 
and 



?Mqj 



(3.9) 



A rr ,(^,q) = ^^J^?7r(p)?7r'(p)x s ( 

"VP 

x\s (-iu>i> - iui, -p- 



(3.10) 



In Eq. (|3.10[) . two Xs's appear because of the dual- 
spin operator. It should be noted that 2Q's are 
equivalent to (±n/2a,0) and (0,±7r/2a): 2Q - G = 
(±7r/2a,0) and (0,±7r/2a), with G = (±27r/a,0) and 
(0,±27r/a) being reciprocal lattice vectors in 2DBZ. 
Then, Cu-0 bond stretching modes around (±7r/2a, 0) 
and (0,±7r/2a) in 2DBZ can be soft provided that 
AF fluctuations around Q = (±37r/4a, ±7r/a) and 
(±7r/a, ±37r/4a) in 2DBZ are developed. 

Since Cu-0 bond stretching modes around 2Q are con- 
sidered, the vibrations of Cu ions are ignored, that is, it 
is assumed that 



\CdVd,\^M p /M d = 0, 



(3.11) 



and 



|C> p , Aq | = Cp eV/A, (3.12) 

where c p is a dimensionless constant and it is likely 36 

c p = 0{l). (3.13) 

Since the contribution from small p is large in the sum- 
mation over p in Eq. (|3.10p . only the contribution from 
the r = s channel is considered. Then, it follows that 



S A ph) K q + J 0,q) 



with 
and 



2M p w Aq 4 2 



3 r AF [ Xs (o,Q)k 2 ] 2 \c p v p . m \ 



(3.14) 



(3.15) 



(3.8) 



s(w Aq + jO,q) = r ?s (q)r ?s (iq) H 2 ■ (3.16) 

Taf [Xs{Q)k. 2 \ 

It should be noted that 3(iu>i, q) is defined as a dimen- 
sionless quantity. The effective transfer integral between 
nearest neighbors for the Gutzwiller quasiparticles is 

t*=t/<j> 7 . (3.17) 

According to Eq. (jl.l| . a plausible number for t* is 

~ W*/8 ~ 40-50 meV. (3.18) 



7 



(a) (iv) (b) (iv) 




k 2 (q x +2n/a)/(n/a) 



FIG. 1: (a) Re[s(cjAq + iO, 2Qo)l as a function of k 2 , with Qo = (— 37r/4a, 7r/a), and (b) Re[H(ojA q + iO, q)] as a function of 
(— 2iv/a < q x < —ir/a) for q y = 2-k j a. For the anisotropy factor, (i) 5 = 1, (ii) 8 — 10~ 1//2 , (iii) 8 = IO" 1 , and (iv) 8 — 10~ 3 . 
In each figure, solid, dotted, dashed, and dashed chain lines are for o^/Faf = 0.2, 0.4, 0.8, and 1.6, respectively. 



According to a microscopic calculation for the spin 
susceptibility, it follows that Taf/K*! = 0(1) and 
Xs(0, Q)k 2 |£*| = 0(1). It is assumed, for the sake of sim- 
plicity, that the energy of Cu-0 bond stretching modes 
is constant and is as large as 

cj Aq = 50 meV. (3.19) 

Then, Aq defined by Eq. (|3 . 1 5|) is approximately given 
by 

Aq ~ 10xc^ M0,Q)K 2 |i*|] 2 meV 

~ I0c 2 p meV. (3.20) 

In this paper, T = K is assumed in the lov sum of 
Eq. (I3TU1) . 

The softening around one of 2Q's or 2Qo, with Qo = 
(— 37r/4a, n/a) in 2DBZ is considered; 2Q is equivalent 
to (jr/2a, 0). Figure [1] shows the dependence of S(^A q + 
i0, q) on k, 5, Taf, and q; Fig. Q|a) and Fig. QJb) show 
Re[S(cjA q + iO, 2Q )] as a function of k 2 and Re[S(aA q + 
iO, q)] as a function of q x , respectively, for several sets of 
S and WAq/TAF- According to Fig. QJb), S(wA q + iO, q) 

has a maximum, that is, Re [E^ ph ^ (w^q + *0, q)] has a 
minimum around 2Qo as a function of q. According to 
Eq. 1(530)1 . Fig. (Ha), and Fig.DJb), it is likely that the 
softening at 2Qo is as large as —(10-20) meV for k 2 «1 
and 5 <C 1. It should be noted that the softening can only 
be large provided that k 2 <C 1 and S <§; 1, as is implied 
by Mermin and Wagner's argument^ integrated effects 
on the softening are never divergent even in the limit of 
k — ► and 8 — * 0. 

It is definite that k 2 <C 1 in the critical region of anti- 
ferromagnetic ordering or spin density wave (SDW), and 
it is certain that the anisotropy is as large as S < 10~ 3 
in cuprate oxide superconductors. Then, the second- 
harmonic effect of AF spin fluctuations can explain the 
observed softening-^ as large as —(10-20) meV around 
2Q or 2Q. 



Since the softening is small when k 2 is large or AF spin 
fluctuations are not developed, it must be small in over- 
doped cuprate oxide superconductors, whose doping con- 
centrations are larger than those of optimal-doped ones. 
When AF spin fluctuations are developed similarly or 
differently between (±37r/4a, ±n/a) and (±7r/a, ±37r/4a) 
because of the anisotropy of the Fermi surface within 
2DBZ, the softening must also occurs similarly or differ- 
ently between the x and y axes or between (±7r/2a, 0) and 
(0, ±7r/2a). These two predictions are consistent with 
observations. 7,8 



IV. STRIPES AND CHECKER BOARDS 

Since the 4a-period and 4a x 4a-period correspond to 
2Q, with Q = (±37r/4a,±7r/a) and [±n/a, ±3?r/4a), a 
plausible scenario for the stripes and checker boards is 
that the complete softening is followed by the stabiliza- 
tion of CDW with 2Q. In general, the 2Q component of 
the density of states, |02q(e)> as a function of e is com- 
posed of symmetric and asymmetric components with re- 
spect to the chemical potential or e — 0. The asymmetric 
component is large when CDW with 2Q is stabilized as 
a fundamental 2Q effect. According to an experiment; 37 
the symmetric component is larger than the asymmetric 
one, which contradicts the scenario of CDW even if the 
softening of the 2Q modes is large and the 2Q fluctua- 
tions are well developed. On the other hand, the symmet- 
ric component is large when the 2Q modulation is due 
to a simple second-harmonic effect of an ordered SDW 
with Q. 51 The second-harmonic effect of the SDW can 
explain the observed almost symmetric P2q(e)- When 
stripes and checker-boards are really static orders, stripes 
must be due to single- Q SDW and checker-boards must 
be due to doublc-Q SDW. It is predicted that magne- 
tizations of the two waves must be orthogonal to each 
other in double-Q SDW . 52 ' 53 It is interesting to examine 
whether the prediction actually holds in cuprate oxide 
superconductors. 
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The appearance or stabilization of SDW is a transi- 
tion rather than a crossover. However, no specific heat 
anomaly has been reported so far except for the anomaly 
due to superconductivity. The absence of any specific 
heat anomaly implies that, even if SDW is stabilized, 
SDW is never a homogeneous phase but is an inhomo- 
geneous phase, which is composed of many domains. If 
the transition temperatures of SDW can be different in 
different domains, no significant specific anomaly can be 
observed. It is plausible that SDW is an disorder-induced 
SDW41 

On the other hand, it is proposed^ that a stripe or 
a checker-board at rather high temperatures must be an 
exotic ordered state, that is, a fluctuating state in a quan- 
tum disordered phase. It should be examined whether 
it is actually such an exotic state. Another possibil- 
ity is that it is a rather normal low-energy fluctuating 
state, whose energy scale is as small as that of the soft 
phonons. The other possibility is the disorder-induced 
SDW^i which can behave as a fluctuating state because 
it is inhomogeneous. 

V. ATTRACTIVE INTERACTION 

Although the electron-phonon interaction plays no or 
only a minor role in the formation of d7-wave Cooper 
pairs, as is discussed in the Introduction, isotope shifts 
of T c can arise from the depression of superconductivity 
by the 2Q fluctuations, whose development depends on 
the mass of O ions. 

In Kondo-lattice theory, cuprate oxide superconduc- 
tors can be relevantly treated as one of the typical Kondo 
lattices. According to Eq. (|2.32[) . which is one of the most 
crucial results of Kondo-lattice theory, two mechanisms 
of attractive interactions, the spin- fluctuation mechanism 
and the exchange-interaction mechanism, are essentially 
the same as each other. However, the attractive inter- 
action mediated by low-energy spin fluctuations such as 
those observed by low-energy neutron inelastic scatter- 
ings or those described by the phenomcnological spin 
susceptibility (|2.24[) is physically different from the at- 
tractive interaction due to the superexchange interaction, 
which arises from the exchange of a pair excitation of 



electrons in spin channels across the Hubbard gap be- 
cause the energy scales of spin fluctuations or spin ex- 
citations are totally different from each other in the two 
physical processes. The main part of the attractive inter- 
action in cuprate oxide superconductors must be the su- 
perexchange interaction rather than the interaction me- 
diated by low-energy AF spin fluctuations. Since the 
superexchange interaction is as strong as J = —(0.10- 
0.15) eVj22 observed high T c can be easily reproduced, 
as is discussed in the Introduction. 



VI. CONCLUSION 



In cuprate oxides superconductors, the electron- 
phonon interaction arising from the modulation of the 
superexchange interaction by lattice vibrations is strong 
enough to cause the softening of not only the half- 
breathing modes around (±7r/a,0) and (0, ±7r/a) in 
the two-dimensional Brillouin zone but also Cu-0 bond 
stretching modes around (±7r/2a,0) and (0,±7r/2a). Al- 
though the softening of the bond stretching modes is 
responsible for stripe and checker-board fluctuations in 
charge channels, the stabilization of a charge density 
wave state following the complete softening of the bond 
stretching modes can never be any relevant scenario 
for ordered stripe and checker-board states. The or- 
dered stripe or checker-board state must be simply a 
single-Q or double-Q spin density wave state, whose Q's 
are (±37r/4a, ±7r/a) and (±w/a, ±37r/4a). The strong 
electron-phonon interaction can play no or only a minor 
role in the binding of c?7-wave Cooper pairs in cuprate 
oxide superconductors, because the attractive interac- 
tion arising from the virtual exchange of a phonon is 
never strong between quasi-particles on nearest-neighbor 
Cu ions on a Cu02 plane. However, isotope shifts of T c 
can arise from the depression of superconductivity by the 
stripe or checker-board fluctuations. Since the superex- 
change interaction is as strong as J = —(0.10-0.15) eV 
between nearest-neighbor Cu ions, the superexchange in- 
teraction must be mainly responsible for the binding of 
the ci7-wave Cooper pairs in cuprate oxide superconduc- 
tors. 
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